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In  this  paper,  a  time  -  dependent  model  of  prophylactic  vaccination  HBV 
is  considered.  This  epidemic  model  is  being  investigated  for  various 
infectious  diseases  classes,  counsel  and  uncounsel.  Mathematical 
analyzes  are  performed  to  determine  the  positivity  and  we  applied  an 
optimal  control  strategies  in  the  form  of  vaccination  to  minimize  or 
eradicate  transmission  from  mother  to  child.  The  study  concluded  that 
prophylactic  vaccination  is  not  an  efficient  way  to  curb  the  epidemic. 


1.  Introduction 

Hepatitis  is  an  inflammation  of  the  liver  caused  by  viruses,  bacterial  infections,  or  constant 
vulnerability  to  alcohol  intake,  drugs,  or  toxic  chemicals,  such  as  those  obtained  in  aerosol  sprays 
and  paint  thinners  [1].  Inflammation  is  the  painful,  red  growth  that  results  when  tissues  of  the  body 
become  wounded  or  infected.  Inflammation  can  cause  organs  not  to  work  properly.  Hepatitis  can 
also  result  from  an  autoimmune  disorder,  in  which  the  body  wrongly  sends  disease-fighting  cells 
to  attack  its  healthy  tissue,  in  this  case,  the  liver  [2].  The  liver  is  located  in  the  upper  right-hand 
side  of  the  abdomen,  mostly  behind  the  rib  cage.  The  liver  of  an  adult  typically  weighs  close  to 
three  pounds  [3].  No  matter  its  cause,  hepatitis  decreases  the  liver’s  ability  to  make  life¬ 
preserving  functions,  including  filtering  fatal  infectious  agents  from  the  blood,  storing  blood  sugar 
and  converting  it  to  usable  energy  forms,  and  producing  many  proteins  necessary  for  life  [4]. 
World  Health  Organization  [5]  has  described  Hepatitis  B  as  a  potentially  life-threatening  liver 
infection  of  global  importance.  An  approximated  240  million  people  worldwide  are  chronically 
affected  with  this  virus,  and  more  than  780,000  people  die  each  year  from  the  complexities 
associated  with  Hepatitis  B.  The  infection  occurs  in  two  phases.  The  first  is  an  acute  phase,  which, 
in  most  people,  doesn’t  result  in  any  ill-health  symptoms.  However,  a  few  people  may  experience 
an  acute  illness  lasting  several  weeks  with  symptoms  including  jaundice,  extreme  fatigue,  dark 
urine,  nausea,  vomiting  and  abdominal  pain.  Fewer  still  can  develop  acute  liver  failure,  which  can 
lead  to  death.  The  second  phase  is  a  chronic  carrier  state  in  which  the  virus  remains  in  the  liver 
and  may  ultimately  cause  chronic  liver  insufficiency,  liver  cirrhosis  or  primary  liver  cell  cancer. 
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These  carriers  show  no  overt  signs  of  ill-health,  and  they  are  usually  the  source  of  all  new 
infections. 

2.  Mathematical  Model  Formulation 

We  propose  a  mathematical  model  to  understand  the  transmission  dynamics  and  prevalence  of 
HBV.  The  model  is  constructed  based  on  the  characteristics  of  HBV  transmission  and  the  model 
[6].  Note  that  they  proposed  a  six  compartmental  model  for  the  propagation  of  the  epidemic 
model.  It  is  observed  that  the  immunity  after  recovery  is  not  lifetime,  while  the  vaccination  might 
wane  after  some  time  [7].  Figure  1  show  the  schematic  representation  of  interactions  of  HBV 
transmission.  Therefore,  we  divide  the  host  population  into  seven  epidemiological  groups:  the 
proportion  susceptible  to  infection  S;  those  latently  infected  E;  uneducated  acute  infections  E; 
educated  acute  infections  le;  carriers  C;  recovered  R;  and  Vaccinated  V. 


Figure  1:  Schematic  representation  of  interactions  of  HBV  transmission 
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In  this  model,  the  Holling-type  II  incidence  function  given  by  g  (/)  = - . 

\^- col 

The  choice  of  this  incidence  function  is  necessary  due  to  the  preventive  measure  (and  behavioral 
changes)  taken  by  the  susceptible  individuals  in  response  to  the  severity  of  the  disease. 

The  parameters  used  are  defined  in  Table  1 . 


Table  1;  the  parameters  of  the  model  are  defined  in  the  table  below 


Parameter 

Interpretation 

Value 

Reference 

Birth  rate 

0.0121 

[8] 

Mq 

Natural  mortality  rate 

0.00693 

[8] 

Mx 

HBV  related  mortality  rate 

0.007 

[9] 

P 

Transmission  coefficient  or  effective  contact  rate 

0.00131 

[8] 

CO 

Proportion  of  births  without  successful  vaccination 

0-100  percent 

[61 

q 

Average  probability  an  individual  fails  to  clear  an  acute  infection  and 
develops  to  carrier  state 

0.885 

[10] 

Yo 

Rate  of  movement  from  uncounsel  acute  to  carrier 

4  per  year 

[7] 

Yi 

Rate  of  movement  from  counsel  acute  to  carrier 

2  per  year 

[7] 

Yi 

Rate  of  movement  from  carrier  to  immune 

0.025  per  year 

[7] 

Y3 

Vaccination  rate  of  susceptible 

0-100  percent 

[7] 

Y4 

Rate  of  movement  from  uncounsel  acute  to  counsel  acute 

Hypothetical 

cr 

Rate  of  movement  from  latent  to  acute 

6  per  year 

[7] 

S 

Reduced  transmission  rate 

0.16 

[7] 

¥ 

Rate  of  waning  of  vaccine-induced  immunity 

0.04 

[11] 

V 

Proportion  of  prenatally  infected  (carrier  mothers) 

0.11 

[7] 

The  total  population  N  (t)  can  be  obtained  from 


(v(o=s(0+£(0+4(0+^.(0+c(t)+/?(t)+v(t) 


(9) 

(10) 


Here,  it  is  important  to  note  that  in  the  absence  of  the  disease  Nit)  — >  —  Moreover,  under  the 

/^o 

dynamics  described  by  the  above  systems  of  equations,  the  region. 
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Q=  Y  =  {S,E,I^J^,C,R,  V,N)^^l\S>Q,E>0,  /„  >0,7,  >0,C>0,y  >0,iV<^|  (ll) 

is  said  to  be  positively  invariant.  Hence  the  system  is  both  mathematically  and  epidemiologically 
well-posed.  Therefore,  for  initial  starting  point  the  trajectory  lies  in  Q.  Thus  we  restrict 

our  analysis  to  the  region  Q.  (Where  the  models  make  biological  sense) 

Theorem  1  All  the  solution  of  the  equations  (1-7)  are  positive  for  all  time  t  >  0  provided  the 
initial  condition  are  positive. 

Proof:  Let  {(5(0),  £(0),  7„(0),7,(0),C(0),/?(0),L(0))>0}eiR’ 


—  =  fio){\-vC)  +  \f/V -{X  +  fi^+el3C  +  r,)S{t) 

^  — (/i  +  /Uq  +  spc  +  )  5(r) 

This  implies,  5'(r)  > -{^X  +  fdf^  +  sfiC  +  y^^Sit),  integrating  we  have 


5(r)>5(0)e 


-j{A+£fiC)dt-{fi„+r3)t 


>0 


(12) 

(13) 

(14) 


Similarly,  it  can  be  shown  that  5  >  0,  £■  >  0,  /„  >  0,  /^  >  0,  C  >  0,  L  >  0,  i?  >  0, ,  for  all  time  t  >  0. 
Hence  all  solutions  of  the  HBV  model  remain  positive  for  all  non-negative  initial  conditions. 


3.  Optimal  Model 

dS 

—  =  jua)(l  -vC  +  vu^C)  +  y/V  -  (//q  +  1  +  £)5C  +  )  5 

dt 
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dR 
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(15) 

(16) 

(17) 

(18) 

(19) 

(20) 


The  problem  is  to  minimize  or  eradicate  completely  mother  to  child  transmission  and  the  cost 
control  implication  on  the  vaccination  strategies  employed.  The  objective  function  is  given  as; 


7(^1,112)  —  J*  Aj5  (r )  +  A2C (r) + 


-u 


(r)  dt 


(21) 
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Where,  the  rate  of  professional  effort  that  the  medical  expatriate  put  in  place  in  prevention  or 
eradicate  of  mother  to  child  transmission  (mj)  and  cost  of  controlling  the  vaccination  {u^).  The 

total  cost  includes  not  only  the  consumption  for  every  individual,  but  also  the  cost  of  organization, 
management,  and  cooperation,  etc.  Hence,  the  coefficients,  Ai,  A2,  Bi,  and  B2,  are  balancing 
cost  factors  due  to  scales  and  the  importance  of  the  four  parts  of  the  objective  function  above. 
Using  similar  control  variables  and  objective  function  (as  in  above)  we  find  that  the  Hamiltonian 
takes  the  form; 

H  =  A,S{t)  +  A^C{t)  +  ^ul{t)  +  ^ul[t)  +  ^A.J^ 


Where, 

Zo=Ao+c^>  Zi  =//o  +  ro+r4’  h=M(i+rv  Z^=-M(OV  +  jUo+Mi+r2^  Z^=jUo+¥- 
B  B 

H  =  A2C J — ~  1^2  (t) “h dj  (^1  ~ vC  +  vu^C ^ -1- y/V  —  ju^  +  A. -h sfiC -1- Y'^2 ) J 

-1-/I2  ^[X  +  s/3C)S-ZqE~^+X2\p'E-ZJ^^+X^\yJ^  [//®v(l-Mj)C-i-^(yo4 

+  X^^fd{l-co)  +  y2U2S -Z^~^  (23) 


/Ij  —  ~Aj  -|-  /Ij  ^/-Iq  -h  A-h  Sf^C  -f  ^3^2  ]~  £f5C^A2~ 

(24) 

=  Z^A^  ~ 

(25) 

/?3  —  x2 

(1  +  «A) 

(26) 

(l  +  «2/J 

(27) 

A^'  =  -A2  +  //®v(l-Mj)(/l5  - Aj)  +  jU(i)vAj  -£/3S(Aj  )  +  Z3/l5 

(28) 

(29) 

With  transversality  conditions  A^  )  =  0  (i  =  1, ...,  6,  is  the  end  time).  Using  the 
Hamiltonian,  we  obtain  the  optimality  conditions. 

(A.  -A.)//covC 
u,  =  - - - - 

(30) 

1/t^ 

B2 

(31) 
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4.  Results  and  Discussion 
4.1.  Numerical  Solution 

To  obtain  further  insight  into  the  influence  of  the  control  parameters  on  the  behavior  of  the 
systems  under  consideration,  thus,  we  consider  the  following  Figures  plotted  with  the  aid  of  Maple 
18  Software  [12], 


- Carrier 

acute  uncounsel  infection 

Susceptible 

Vaccination 

- Carrier 

acute  uncounsel  infection 

Susceptible 

Vaccination 

Figure  2:  Graph  of  Carrier,  Acute  uncounsel 
infectious.  Susceptible  and  Vaccination  at 

Ui=l, 

initial  stage  (without  control) 


Figure  3:  Graph  of  Carrier,  Acute  uncounsel 
infectious.  Susceptible  and  Vaccination  at 

U2=0.9 


'  Carrier 
Susceptible 


acute  uncounsel  infection 
■  Vaccination 


■  Carrier 
Susceptible  • 


■  acute  uncounsel  infection 

■  Vaccination 


Figure  4:  Graph  of  Carrier,  Acute  uncounsel 
infectious.  Susceptible  and  Vaccination  at 
ui=0.1. 
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Figure  5:  Graph  of  Carrier,  Acute  uncounsel 
infectious.  Susceptible  and  Vaccination  at 
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Ui=0.1,  U2=0 


Carrier  atul=1.0  andii2=0.9  (b&<^ 

. Carrier  at  ul  =0.1  and  u2=0.9 

- Acute  uncounsd  infection  at  ul=l  .0  and  u2=0-9 

. Acute  uncounsd  infection  at  ul=0.1  and  u2=0.9 

Figure  6:  Graph  showing  the  sensitivity 
in  Carrier  and  Acute  uncounsel  infectious  in 
figure  2  and  4 


- Carrier  at  ul=0.7,  u2=0.3 - Acute  uncounsd  infection  at  ul=0.7,  u2=0.3 

- Carrier  at  ul=r  u2=0  Acute  uncounsd  infection  at  ul=l,  u2=0 

- Carrier  atul=0.4  u2=0.6 - Acute  uncounsd  infection  atul=0.4,  u2=0.6 

. Carrier  at  ul=0.L  u2=0.9 . Acute  uncounsd  infection  at  ul  =0.1,  u2=0.9 

Figure  8:  Graph  showing  the  sensitivity 
in  Carrier  and  Acute  uncounsel  infectious 
for  different  values  of  control 


U2=0.9 


- Carrier  at  ul  =1 .0  and  u2=0  (a&c) 

- Acute  uncounsd  infection  at  ul=l  .0  and  u2=0 

. Carrier  at  ul  =0.1  and  u2=0 

Acute  uncounsd  infection  at  ul=0.1  and  u2=0 

Figure  7:  Graph  showing  the  sensitivity 
in  Carrier  and  Acute  uncounsel  infectious  in 
figure  3  and  5 


- Ul=<l,  U2=0  - U1=0,  U2=0.33 

- U1=0,  U2=0.66 - U1=0.  U2=0.99 


Figure  9:  The  plot  shows  the  change  in 
Susceptible  populations  at  different  values 
of  control 
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- U1=0,  U2=0  - U1=0,  U2=0.33 

- U1=0,  U2=0.66 - U1=0,  U2=0.99 


Figure  10:  The  plot  shows  the  change  in  Vaccination 
populations  at  different  values  of  control 


To  understand  the  relationship  between  the  two  control  parameters  (the  rate  of  prevention  or 
eradicates  of  mother  to  child  transmission  ui  and  the  cost  of  controlling  the  vaccination  U2)  in  the 
model,  we  solved  the  optimal  control  problem  when  ui  =  0  and  U2  =  0  (no  control)  (Figure  2),  by 
setting  ui  =  1  and  U2  =  0.9  (Figure  3)  and  looking  at  the  resulting  scale  for  ui  =  0.1  and  U2  =  0 
(Figure  4)  and  also  considering  when  ui  =  0.1  and  U2  =  0.9  (Figure  5).  The  four  scenarios  show 
two  major  significant  difference.  Figure  2  and  4,  Figure  3  and  5,  respectively,  show  small 
difference  as  portrayed  by  Figure  6  and  7.  We  own  this  change  to  the  control  parameter  ui,  while 
the  major  change  (Figure  8)  recorded  is  due  to  the  cost  effectiveness  of  the  vaccination  strategies 
U2.  Moreover,  we  also  show  the  effect  of  the  control  parameters  on  both  the  Susceptible  and 
Vaccination  population  in  Figure  9  and  Figure  10  respectively.  This  delivers  significantly  better 
results  due  to  the  split  in  the  acute  infectious  class  into  acute  counsel  and  acute  uncounsel 
infectious  classes.  However,  when  comparing  our  results  to  those  of  older  studies  [6],  it  must  be 
pointed  out  that  the  acute  class  differs  and  the  optimal  strategies  in  transmission  and  vaccination 
are  better  approach  in  assessing  how  well  the  model  behaves. 

5.  Conclusion 

In  this  paper,  we  proposed  seven  compartmental  models  of  Hepatitis  B  virus  infection  with  two 
controls;  the  prevention  or  eradicates  of  mother  to  child  transmission  (ui)  and  cost  of  controlling 
the  vaccination  strategies  employed  (U2).  The  positivity  of  the  model  is  obtained.  Thus,  we 
formulate  an  optimal  control  problem  and  solved  it  using  the  Pontryagin’s  maximum  principle. 
The  results  established  shows  that  effort  put  in  place  in  subduing  or  eliminating  the  birth  of  carrier 
by  prenatally  infected  carrier  mother  is  not  a  sufficient  intervention  to  employ.  The  study 
concluded  that  prophylactic  vaccination  is  not  a  sufficient  way  to  curb  the  epidemic. 
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